Compiling gate networks on an Ising quantum computer 
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Here we describe a simple mechanical procedure for compiling a quantum gate network into the 
natural gates (pulses and delays) for an Ising quantum computer. The aim is not necessarily to 
generate the most efficient pulse sequence, but rather to develop an efficient compilation algorithm 
that can be easily implemented in large spin systems. The key observation is that it is not always 
necessary to refocus all the undesired couplings in a spin system. Instead the coupling evolution can 
simply be tracked and then corrected at some later time. Although described within the language 
of NMR the algorithm is applicable to any design of quantum computer based on Ising couplings. 

PACS numbers: 03.67.Lx 



Quantum computers [lj have generated considerable 
interest in recent years due to their apparent ability to 
perform computations that are intractable on any clas- 
sical computer. Although the construction of a general 
purpose quantum computer capable of solving real prob- 
lems remains a challenge, preliminary results have been 
demonstrated in several systems 0, most notably nu- 
clear magnetic resonance (NMR) 0,0, |H @ and trapped 
ions 0j la, 0> lldj . When implementing simple quantum 
algorithms on devices with a small number of qubits, 
it is perfectly possible to design an implementation "by 
hand" . With larger systems, however, this approach be- 
comes impractical, and it is desirable to automate the 
process. Here we describe a simple procedure that allows 
a sequence of abstract quantum logic gates (the most 
common device- independent description of a quantum al- 
gorithm) to be compiled into a sequence of pulses and 
delays, which are the natural gates on an Ising quantum 
computer. 

By an Ising quantum computer we mean a system with 
a background Hamiltonian 



(1) 



where, following NMR notation ^lj, I z = a z /2 is the 
angular momentum operator in the z direction, is the 
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precession frequency of the ith qubit, and Jjj is the Ising 
coupling strength between pairs of qubits; here we as- 
sume that all the n(n — l)/2 Ising couplings in an n 
qubit system are of significant size. This Hamiltonian 
arises naturally in NMR systems, where the qubits cor- 
respond to spin-half particles and the scalar J coupling 
takes the Ising form within the weak coupling limit, but 
also occurs more generally. It can be simplified by work- 
ing in a multiply rotating frame, where the frame used to 
describe each qubit rotates at the Larmor precession fre- 
quency Vi, so that the effective Hamiltonian contains only 
the coupling terms. We also assume that the system can 
be controlled by the application of arbitrary single-qubit 
gates to any qubit or group of qubits. In an NMR system 
these single-qubit gates would be implemented using res- 
onant radiofrequency pulses, with individual qubits being 
distinguished by their unique Larmor frequencies. We re- 
fer to these single-qubit gates as pulses whether or not 
they are implemented in this way. 

It is well known that any quantum algorithm can be 
implemented using only single-qubit and two-qubit gates, 
such as CNOT (controlled-NOT) gates . Further- 
more, the CNOT gate can itself be decomposed in terms of 
(single-qubit) Hadamard gates and a controlled-cr z gate, 
which converts |11) to — 111) while leaving other states 
unchanged. Finally, a controlled-cr z is essentially equiva- 
lent to a 90° evolution under an Ising coupling, differing 
only by single-qubit z rotations |5|. These are most con- 
veniently handled not by rotating the qubit, but rather 
by rotating its reference frame, a method usually known 
as abstract reference frames [l4j . Similar techniques can 
be used to implement more general controlled phase-shift 
gates , which are useful for building such gates as the con- 
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FIG. 1: An example of rescaling couplings between four 
qubits. The solid blocks indicate NOT gates (180° pulses). 



trolled square-root-of-NOT |5|. Thus arbitrary quantum 
circuits can be constructed out of pulses and controlled 
evolutions under the Ising Hamiltonian. Here we simply 
assume that a circuit comprising only single-qubit gates 
and controlled phase-shift gates is available. 

In order to achieve the desired evolutions it is necessary 
to sculpt the Hamiltonian into a more suitable form, in 
which only the desired couplings of the required strength 
are present. This is easily achieved by using spin echoes. 
Spin echoes are widely used in conventional NMR experi- 
ments to refocus undesirable single-spin interactions, 
but they can also be used to effectively rescale Ising cou- 
pling terms within a Hamiltonian. 

The basic procedure, which is based on the concept of 
an average Hamiltonian, is easily understood and is ex- 
emplified by Fig. ^ The coupling between qubits and 
1 evolves for the entire period, To, and so the average 
coupling strength experienced by these qubits is equal 
to the underlying coupling strength, J i. The coupling 
from qubit to qubit 2 is partially refocused by applying 
a pair of not gates (180° pulses) to one of the two qubits 
involved. The effect of these not gates is to reverse the 
coupling evolution Hi for a period of length e2, and 
so it appears that the coupling has evolved only for a time 
T2 = To — 2e2. Alternatively (and equivalently) the situ- 
ation can be described as if the coupling had evolved for 
the whole time To but with a reduced coupling strength 
Jq 2 = (T2/T0) Jo2- In the same way the coupling between 
qubits and 3 has been rescaled by a factor of t^/tq, 
with t 3 = t — 2e3. The effective coupling strengths be- 
tween other pairs of qubits can be worked out in a similar 
way. Clearly J' l2 = (t 2 /t ) J12 and J' 13 = (t 3 /t ) J13, by 
analogy with the couplings to qubit 0. The last prob- 
lem is to calculate the net evolution for the coupling be- 
tween qubits 2 and 3, where NOT gates are applied to both 
qubits, and a little thought shows that the scaling factor 
is T23/V0 with a net evolution time T23 = To — 2(63 — £2). 

The standard approach to date has been to use refo- 
cusing networks to sculpt the underlying Ising Hamilto- 
nian into some ideal form, usually retaining only a small 
number of coupling interactions. In p articular, efficient 
refocusing schemes are known [la . Il6| to retain any one 
coupling in an extended network. There are, however, 



some limitations on the average Hamiltonians that can 
be achieved. For example, consider a system of three 
coupled spins, with a Hamiltonian containing three cou- 
pling terms: it is possible to keep any one of these terms 
while refocusing the other two, but it is not possible to 
keep two terms while refocusing the third. If such an 
average Hamiltonian is desired it is usually necessary to 
implement the two Hamiltonians sequentially. Designing 
evolution networks in large systems can become a com- 
plex business. 

Our compilation algorithm provides an alternative ap- 
proach for designing pulse sequences on Ising quantum 
computers with large numbers of spins. The aim is not 
to generate the most efficient pulse sequence, but rather 
to develop a simple algorithm that can be easily imple- 
mented in large spin systems. The key observation is that 
it is not, in fact, necessary to implement the exact Ising 
coupling evolutions shown in the network. In particu- 
lar, it is not always necessary to refocus all the undesired 
Ising couplings in a spin system. Instead the coupling- 
evolution can simply be tracked, and their values cor- 
rected at some appropriate time. This is because single- 
qubit gates applied to one qubit commute with couplings 
that do not involve this qubit. Thus it is only neces- 
sary to achieve the "correct" evolution for those qubits to 
which single-qubit gates are applied. Furthermore, since 
single-qubit gates applied to different qubits all commute 
with one another, two or more "simultaneous" single- 
qubit gates can in fact be applied sequentially in any 
order. 

The first stage is to redraw the network so that single- 
qubit gates are applied sequentially, rather than in par- 
allel. Simultaneous gates can be applied in any order; 
the simplest approach is to apply them from top to bot- 
tom. We refer to the qubit to which a single-qubit gate 
is applied as the target qubit, and any other qubits as 
control qubits. It is not necessary that evolution periods 
generate any particular coupling evolution between the 
control qubits, but these additional evolutions must be 
tracked through the circuit. 

Couplings are tracked by recording the net evolution 
angle (modulo 360°) generated for each of the n(n — 1) /2 
couplings at all significant points in the network, where 
the change in the evolution angle is given by 

<5% = TT^.T, (2) 

and J[j is the average coupling strength between qubits 
i and j during a time period t. It is essential that all the 
coupling angles involving a target qubit have the correct 
value as defined by the gate network before applying a 
single-qubit gate to the target. This is easily achieved by 
applying not gates to the control spins, thereby changing 
the average coupling strengths, so that all the couplings 
reach the desired angle in the same evolution time. When 
the single-qubit gate is applied to the target qubit, all 
coupling angles to this qubit are reset to zero. Couplings 
between control qubits are simply tracked throughout the 
process. 
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FIG. 2: An example quantum circuit, based on the algorith- 
mic benchmark of Knill et al. [T^|. White squares indicate 
90°, pulses (pseudo-Hadamard gates) while black circles con- 
nected by control lines are Ising coupling gates with a target 
angle of 90° between the coupled qubits. 



TABLE I: Coupling evolution angles (measured in degrees, 
and rounded to the nearest degree) between every pair of 
qubits at various points in the initial part of the circuit shown 
in Fig. |21 Coupling angles indicated in boldface are coupling 
angles to target qubits that must be set correctly; other angles 
are simply being tracked as part of the algorithm. 
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FIG. 3: An example Ising quantum computer, described by 
indicating the strengths of the Ising couplings (measured in 
Hz) between all pairs of qubits; these couplings were chosen 
arbitrarily. 



This procedure can be clarified by considering a con- 
crete example, and here we consider the implementation 
of the quantum circuit shown in Fig.[21on the Ising quan- 
tum computer depicted in Fig. [21 The circuit is based on 
the algorithmic benchmark of Knill et al. but the 

coupling strengths shown in the Ising computer were cho- 
sen at random. The key values in an Ising computer are 
the coupling angles achieved between pairs of qubits, and 
these will be considered at the points marked with lower 
case letters in Fig. 

The first part of the circuit is simple: the system be- 
gins with all coupling angles set to zero, and the first 
single-qubit gate has no effect, so that at point (a) all 
coupling angles are 0°. By point (b) it is necessary that 
the coupling angle between qubits and 1 be increased to 
90°, while all other coupling angles to qubit 1 (which is 
the target of the next single-qubit gate) remain at zero; 
all further couplings can simply be tracked. This can 
be achieved by evolving under the coupling Hamiltonian 
for a time l/2Joi = 11905 /is with NOT gates applied to 
qubits 2 and 3 at the middle and end of the coupling pe- 
riods, acting to completely suppress the couplings from 
these qubits to qubit 1. The couplings from the other 
control qubits (2 and 3) to qubit are also suppressed, 
but the coupling between qubits 2 and 3 will evolve for 



the whole coupling time, and thus achieve a coupling an- 
gle of 144°. The next step, at point (c), is to apply a 
single-qubit gate to qubit 1; at this point all coupling 
angles involving this qubit are reset to zero, while leav- 
ing all other angles unchanged. Thus the three coupling- 
angles involving only control qubits still have the values 
#02 = 0°, 6» 03 = 0°, and 8 23 = 144°. These (and further) 
evolutions are summarized in Table [I] 

The first interesting point in the algorithm occurs at 
(d). At this stage the algorithm requires the coupling 
angles to qubit 2 to be set correctly, that is 8q2 = 0°, 
812 = 90°, and 823 — 0°. Allowing for the current val- 
ues of the coupling angles, and recalling that coupling 
angles are defined modulo 360°, the additional evolution 
required is <5 02 = 0°, 5 12 = 90°, and S 23 = 216°. Con- 
sidering these evolutions individually, the coupling times 
required are T02 = /is, t\% — 8621 /is and T23 = 17910 /is. 
Clearly the last coupling requires the longest time, and so 
all the desired coupling angles can be achieved by evolv- 
ing under the full coupling Hamiltonian for 8955 /is, with 
two couplings, J02 and J12, being (respectively) totally 
and partially refocused. Thus not gates should be ap- 
plied to qubit at the middle and end of the coupling 
period, while the first NOT gate is applied to qubit 1 after 
a time T12 + {j2z — i"i2)/2 = 13266 /is, with the second 
not gate applied at the end. The remaining coupling 
angles can then be calculated by noting that J03 is com- 
pletely refocused, J13 evolves for the same time as J12, 
giving an angle of 76°, and J01 evolves for a net time of 
9288 /is, giving an additional coupling angle of 70° . 

At point (e), a single-qubit gate is applied to qubit 
2, and so all coupling angles involving this qubit are re- 
set to zero. The couplings at point (f) are then set using 
the same method as was used for (d). This time the limit 
coupling is J13 which requires a time of T13 = 32200 /is to 
acquire an additional coupling angle of 813 — 284°. The 
coupling J03 is completely refocused by applying not 
gates to qubit at the middle and end of this coupling 
period, while the coupling J23 must be partially refo- 
cussed by applying the first NOT gate to qubit 2 after a 
time of 19832 /is. The remaining couplings are tracked 
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FIG. 4: Sequence of refocusing pulses used to implement the circuit shown in Fig. [5] on the Ising quantum computer shown 
in Fig. |2l White squares indicate 90y pulses as before, while white rectangles indicate NOT gates. These single-qubit gates 
are assumed to be instantaneous. The evolution periods are drawn approximately to scale, and the lengths of some periods 
(measured in fjs) are indicated. 



as usual, and their values are shown in Table H The 
complete sequence of refocusing pulses and single-qubit 
gates is shown in Fig. Note that the evolution times 
depicted are those calculated above, which include the 
effects of rounding errors; exact calculations would give 
slightly different times. 

As discussed above, this algorithm is not necessarily 
intended to produce the most efficient pulse sequence, 
but it is still necessary to check that the simplicity of im- 
plementation is not achieved at an excessive cost either in 
pulse sequence length or in the number of pulses required. 
Clearly the complexity of a pulse sequence will depend 
on details of the circuit being implemented and the range 
of coupling strengths in the experimental system. It is, 
however, simple to consider some limiting situations. 

Any portion of a circuit in a system of n qubits can 
be classified according to the number p of target qubits 
involved, which obviously lies between the extreme val- 
ues of 1 and n. Using our new algorithm this section 
will be divided into at most p separate sections, each of 
which will require 2(n — 2) refocusing pulses; in total the 
implementation will require less than 2xnxp pulses. Us- 
ing the Hadamard matrix approach will require between 
1 and p refocusing periods (depending on exactly which 
couplings have to be controlled), each of which will con- 
tain approximately n 2 refocusing pulses [Tc| . The relative 
efficiency of our new algorithm against the Hadamard al- 
gorithm will depend on the number of refocusing periods 
required by the latter. In a system such as that shown in 
Fig. |3 where all the coupling strengths are different, it 
will not be possible in general to combine different evolu- 
tion periods, and so the Hadamard approach will require 
about p x n 2 pulses. Thus our new algorithm will nor- 
mally use fewer pulses than previous methods. The time 
required to implement a set of gates is another important 
consideration, but once again will vary greatly from net- 
work to network. It is easy to sec that pulse sequences 
produced by our new algorithm will in the worst case take 
p times longer than those produced by current methods, 
although it seems likely that the relative performance will 
be better than this in real situations. 

It is possible to imagine a large number of ways in 



which this algorithm could be extended, resulting in sim- 
pler or shorter pulse sequences. Here we confine ourselves 
to two particularly simple extensions, both of which can 
be easily implemented. Firstly, we note that it is not in 
fact necessary to place our pairs of not gates at the end 
of the coupling period, as shown in Fig. ^ Instead the 
pair of refocusing gates can be placed at any point within 
the coupling period without affecting the net coupling to 
the target spin (the position will, of course, affect the 
net couplings to other control spins). In particular, we 
can choose to place the first NOT gate at the beginning 
of a coupling period, and if this period is preceded by 
one with a NOT gate at the end then two of the four not 
gates can be canceled. 

Secondly, as described so far, it has been necessary to 
set coupling evolution angles to a target spin to the cor- 
rect value modulo 360°; thus qubit pairs can, in principle, 
be required to undergo evolutions of up to 360° to achieve 
the desired angle. This is, however, excessive. We note 
that 

exp[-*7r 2P z I k z ] = expH^i + J*)] (3) 

and so a 180° coupling evolution on a pair of qubits 
is equivalent to a 180° frame rotation applied to both 
qubits. Thus it is never necessary for qubits to undergo 
an evolution through more than 180°, as any additional 
evolution can be achieved through frame rotations. 

Thirdly, we note that, in addition to partly refocusing 
coupling periods, not gates can also be used to negate 
the sign of an Ising coupling. If a pair of not gates 
is applied to a control qubit at the beginning and end 
of an evolution period, then the evolution of the Ising 
coupling will be reversed throughout this period. If NOT 
gates are already applied to the qubit during this period 
to partially refocus a coupling, then two of the four not 
gates can be cancelled, so there will be no overall increase 
in the number of gates. Combined with the previous 
observation, we note that any desired evolution can be 
achieved using a coupling evolution in the range ±90°. 
This allows the implementation time to be reduced by a 
factor of around four, at the cost of a small increase in 
the number of not gates. 
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Finally, we note that it is possible to combine the tra- 
ditional approach and our new algorithm to produce hy- 
brid schemes. For example, it might be useful to compile 
portions of a network independently of any larger com- 
putation, in effect treating them as subroutines. To do 
this, it is necessary to ensure that all the couplings in 
the system have correct evolution angles at the end of 
the subroutine, not just those involving the final target 
qubit. This can be achieved by ending the network with 
a conventional refocusing period. If, however, the whole 
network is compiled together, then this final refocusing 
period will normally not be required. Computations usu- 
ally end with measurements on some of the qubits in the 
computational basis, and these eigenstates commute with 



the Ising Hamiltonian. The results obtained from mea- 
surements of these qubits are, therefore, independent of 
their couplings angles, and it is not necessary to refocus 
couplings to them. 
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